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Fermi surface in KFe 2 As 2 determined via de Haas-van Alphen oscillation 
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We have completely determined the Fermi surface in KFe2As2 via de Haas-van Alphen (dHvA) 
measurements. Fundamental frequencies e, a, f, and j3 are observed in KFe2As2. The first one is 
attributed to a hole cylinder near the X point of the Brillouin zone, while the others to hole cylinders 
at the r point. We also observe magnetic breakdown frequencies between a and ( and suggest a 
plausible explanation for them. The experimental frequencies show deviations from frequencies 
predicted by band structure calculations. Large effective masses up to 19 m e for B || c have been 
found, m e being the free electron mass. The carrier number and Sommerfeld coefficient of the specific 
heat are estimated to be 1.01 - 1.03 holes per formula unit and 82 - 94 mJmol -1 K -2 , respectively, 
which are consistent with the chemical stoichiometry and a direct measure of 93 mJmol -1 K~ 2 [H. 
Fukazawa et al., J. Phys. Soc. Jpn. 80SA, SA118 (2011)]. The Sommerfeld coefficient is about 
9 times enhanced over a band value, suggesting the importance of low-energy spin and/or orbital 
fluctuations, and places KFe2As2 among strongly correlated metals. We have also performed dHvA 
measurements on Bao.o7Ko.93Fe2As2 and have observed the a and /3 frequencies. 



I. INTRODUCTION 

Since the discovery of superconductivity at T c = 
26 K in LaFeAs(0, F) by Kamihara et al., 1 iron- 
pnicitde/selenide high-T c superconductivity has been 
a center of activity in the condensed matter physics 
community. 2 Many different materials have been synthe- 
sized and T c has quickly been raised up to T c ~ 55 K, 3 ~ 6 
The superconducting pairing mechanism is still under 
debate^ spin fluctuations on the one hand and orbital 
fluctuations on the other hand. Basically, the former 
approach predicts an s± gap that changes sign between 
electron and hole Fermi surface (FS) pockets*^ while the 
latter predicts an s ++ gap without sign change. 10 

Bai_ x K :! ;Fe2As2 (Ref. \lu) is one of the most stud- 
ied systems. High-quality single crystals can be grown 
by flux method. The parent material BaFe2As2 or- 
ders antiferromagnetically below the Neel temperature 
Tn = 140 K (Ref. Qjl) and is a moderately correlated 
metal with mass enhancements m* /rriband of 2 - 3 in the 
ground state, m* and rriband being the effective and the 
band mass, respectively. 13,14 As K is substituted for Ba, 
the antiferromagnetism is suppressed and disappears at 
x ~ 0.2, and superconductivity appears^ T c reaches 38 
K at optimal doping x ~ 0.4^i. T c does not disappear till 
X = 1; 15 ~ 17 KFe2As2 is a superconductor with T c = 3.4 

Intriguingly, the superconducting gap structure seems 
to change with x in Bai_ :E K : j;Fe2As2. Near the opti- 
mal doping, a fully gapped s-wave superconductivity has 
been indicated by penetration depth, 19 specific heat, 20,21 



thermal conduct ivity^ and NMR measurements^ On 
the other hand, a nodal gap structure in KFe2As2 has 
been suggested by NMR, specific heat^ 4 - penetration 
depth 2 ^ and thermal conductivity measurements!^ - — A 
laser angle-resolved photoemission spectroscopy (laser 
ARPES) and an NMR study as a function of the compo- 
sition x have found a drastic change in the gap structure 
around x — 0.6 j 29 ' 30 Although it has been under debate 
whether the superconducting state is an s-wave or d-wave 
onej^S the octet-line nodes of the gap recently observed 
in another laser ARPES study^L are compatible with an 
s-wave state with accidental nodes. Theoretical studies 
based on the spin-fluctuations approach suggest a variety 
of gap structures with s- or d-wave symmetries depending 
on band and interaction parameters^ 2 - - — 

In this context, precise determination of the electronic 
structure near the Fermi level Ep in KFe2As2 is highly 
desirable. Previously, we reported de Haas-van Alphen 
(dHvA) measurements on KFe2As2 in a letter^ We ob- 
served one small FS cylinder and two relatively large 
ones. With the aid of band structure calculations, they 
were assigned to a hole cylinder near the X point of the 
Brillouin zone (BZ) and hole cylinders at the T point, 
respectively. We thought that the largest hole cylinder 
at r expected from band structure calculations was not 
observed. Recently, we have observed dHvA oscillations 
in Ba .o7K .93Fe2As2. Interestingly, the largest FS cylin- 
der is clearly observed. Motivated by this observation, 
we have reexamined dHvA oscillations in KFe2As2, and 
have noticed that the largest cylinder is also observed 
in KFe2As2- Thus we have completely determined the 



Fermi surface in KFe2As2 via dHvA measurements. The 
mass enhancement amounts to 9, which is much larger 
than the overall band-width renormalization of 2 esti- 
mated in an early ARPES study 36 and indicates the im- 
portance of low-energy spin and/or orbital fluctuations. 
In addition, we explain magnetic breakdown frequencies 
observed between the a and £ orbits assuming eight mag- 
netic breakdown junctions, the positions of which resem- 
ble those of the octet-line nodes. 



II. EXPERIMENTS, LIFSHITZ-KOSEVICH 
FORMULA, AND YAMAJI MODEL 



factor Rd,t describes the influence of disorder/impurity 
scattering, x* D being the Dingle temperature. We can de- 
termine m* and x* D by fitting Rt iT and Rd,t to T- and 
-B-dependences of experimental oscillation amplitudes at 
constant B and T, respectively. The spin reduction fac- 
tor R s ^ r is due to the interference between oscillations 
from up- and down-spin electrons, and the spin-splitting 
parameter S may be expressed as S = (l/2)g e fffi* with 
an effective g factor g e ff- 

With the present setup and second-harmonic detec- 
tion, amplitudes of oscillations in the detected voltage, 
IV, are related to a r as 



High-quality single crystals of KFe2As2 and 
Ban.n7Kn.Q3Fe2As2 were grown by a self flux method as 
described in Ref. [l8|. Large residual resistivity ratios of 
more than 450 and of ~70 were observed for KFe2As2 
and K .93Bao.o7Fs2As2, respectively^ In the case of 
KFe2As2, nine different crystals were measured and gave 
consistent results. We describe results obtained for the 
most thoroughly measured sample. 

The dHvA measurements were performed in a dilution 
refrigerator and superconducting magnet by using the 
field modulation technique^ The modulation frequency 
and amplitude were mostly / = 67.1 Hz and b = 10.4 mT, 
respectively, and the detection was made at the second 
harmonic (2/). A sample was placed in a balanced pick- 
up coil with its c axis parallel to the coil axis. The field 
direction measured from the c axis is denoted by 9, and 
if necessary a subscript is attached to indicate the field- 
rotation plane. The same setup was used to measure ac 
magnetic susceptibility. 

The dHvA magnetization oscillation M osc due to an 
extremal cyclotron orbit normal to B enclosing the k- 
space area A is given by2£ 
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M osc = 2_^a r sin I — — 
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Here fj,* = m* /m e , m e being the free electron mass. 
The effective mass m* is enhanced over the band mass 
TUband by electron-phonon and electron-electron interac- 
tions. The frequency F is given by F — (h/2Tre)A, and 
(j) r is the phase. Not only the fundamental frequency 
F but also its harmonics (r > 1) appear in M osc . \A"\ 
is the curvature factor: A" — d 2 A/dn 2 , where k is the 
wave number along B. Rp.r is the temperature reduction 
factor, where if is a constant (14.69 T/K). The Dingle 



1 dF 
v r oc(cos9 — ——— sin9)J2{r\)a r , (6) 

b A9 



where Ji is the second-order Bessel function, and A = 
2irFb/B 2 . We distinguish the two amplitudes v r and a r 
when necessary. 

In the case of a purely two-dimensional (2D) electronic 
structure, the FS would be a straight cylinder showing a 
single dHvA frequency F, and Fcos9 would be constant 
irrespective of the field direction 9. With a quadratic in- 
plane dispersion, m* cos 9 would also be constant. Given 
the crystal structure of KFe2As2, each FS cylinder would 
contribute 1.452/1* (mJ mol _1 K -2 ) to the Sommerfeld 
coefficient 7 of the specific heat. 

In reality, there is a c-axis dispersion of the electronic 
band energy, which leads to corrugation of the cylindrical 
FS and produces at least two dHvA frequencies corre- 
sponding to the maximal and minimal FS cross sections. 
Yamaji considered the simplest case where the c-axis dis- 
persion takes a form of cos I c k z , where I c is the interlayer 
distance and equals c/2 in KFe2As2, and derived the an- 
gular dependence of the two frequencies^ 



A p 

F± cos 9 — Fq ± Jq (I c kF tan ( 



(7) 



where Fq and AFb are the average of and the differ- 
ence between the two frequencies at 6 = 0, respectively, 
Jo is the zeroth-order Bessel function, and kp is the 
in-plane Fermi wave number. At magic angles where 
I c kpt&n9 = 7r(n — tp) (n: integer), the two frequen- 
cies coincide, resulting in single-frequency dHvA oscil- 
lations with an enhanced amplitude. It is also known 
that the interlayer magnetoresistance at a constant field 
shows maxima at these field angles as a function of the 
field direction when the field is tilted from the c axis. 
This phenomenon is called angle dependent magnetore- 
sistance oscillation (AMRO) and is used to determine 
the FS in quasi-2D metals. The phase ip is 1/4 for the 
above simplest dispersion but actually varies depending 
on details of the c-axis dispersion. 
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FIG. 1. (color online) Pick-up coil signals detected at the 
modulation frequency as a function of magnetic field for 
KFe2As2 (a) and Bao.o7Ko.93Fe2As2 (b). Upper critical fields 
Bci are indicated. These field sweeps were made on cooling 
the dilution refrigerator. The temperatures when B = B C 2 
were 0.17 and 0.18 K for B || c and B || ab in (a) and 0.22 K 
in (b). 



III. RESULTS 
A. Upper critical field B c2 

Figure [IJa) shows pick-up coil signals detected at the 
modulation frequency / as a function of magnetic held for 
KFe2As2<2^ From these data, we estimate that the upper 
critical field B C 2 in this sample is 1.7 T at T = 0.17 K for 
B || c and 4.9 T at T = 0.18 K for B \\ ab. These values 
are slightly larger than those previously determined for a 
lower-quality sampled We note that the a&-plane critical 
field, which is consistent with Ref. 41, is much lower than 
the orbital critical field (~7 T) estimated from the initial 
slope of B C 2 at T = T c in samples of similar quality to 
the present one ; 42 ' 43 confirming the existence of strong 
spin paramagnetic effects4S 

Figure [IJb) shows pick-up coil signal for B || c in 
Bao.o7Ko.93Fe2As2 as a function of field, from which 
B C 2 in this sample is estimated to be 6.0 T at T = 
0.22 K for B \\ c. A more detailed study on B C 2 in 
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FIG. 2. (color online) Examples of dHvA oscillations for 
KFe2As2 (upper) and Bao.o7Ko.93Fe2As2 (lower). T < 0.1 
K. 
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FIG. 3. (color online) Examples of Fourier transforms in 1/B 
of dHvA oscillations for KFe2As2 (upper three spectra) and 
Bao.o7Ko.93Fe2As2 (lowest). Note that the horizontal axis is 
FcosO. Fundamental frequencies are labeled with Greek let- 
ters. The frequency marked Cu(N) is assigned to the copper 
neck oscillation from copper wire of the pick-up coil. 



Ba .o7K .93Fe 2 As2 is given elsewhere^ 4 - 



B. dHvA oscillations 

Figure [2] shows examples of dHvA oscillations. For 
KFe2As2, the frequency of prominent oscillations for 
B || c is about 2.4 kT both at low fields and at high 
fields, though higher frequencies are also noticeable at 
high fields. The oscillations continue down to below 3 
T.2£ For Bao.o7Ko.93F e 2As2, the frequency of prominent 
oscillations is about 6.7 kT. 
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FIG. 4. (color online) Angular variation of Fourier transforms of dHvA oscillations in KFe2As2 for the (010) plane. The used 
field window is between 10 and 17.65 T. The horizontal axis is Fcos9. The spectra are shifted vertically so that the baseline 
of a spectrum measured at degrees is set at 6. The dashed lines labelled ei t h an d Q.i,k are based on the Yamaji model (see 
text), while those labelled fi,h and j3 are guides to the eye. The right panel shows a frequency region near the ai t h frequencies 
in an expanded scale. The shading indicates a frequency spread of on^ calculated for a ±1° error or distribution of the c axis. 
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FIG. 5. (color online) Same as Fig. 0] but for the (110) plane. 
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FIG. 6. (color online) Angular variation of calculated [(a) and (c)] and experimental dHvA frequencies (b). Note that the 
vertical axis is Fcos8. For the calculated frequencies, both original (a) and adjusted ones (c) are shown (see text). For the 
experimental frequencies (b), the mark sizes are based on the oscillation amplitudes logarithmically. Frequencies assigned to 
fundamentals are shown by filled marks. The solid curves indicating fundamental frequencies in (b) are the same as the dashed 
lines in Figs. [4] and [5] The dotted curves indicate positions of harmonic frequencies calculated from the corresponding solid 
curves. The dashed line labelled Cu-N indicates the copper neck frequency arising from copper wire of the pick-up coil. 



Figure [3] shows representative Fourier transforms of 
dHvA oscillations for the two compounds. Fundamen- 
tal frequencies, e (w , a^ th ), C(l,h), and P, are indi- 
cated. The identification of these fundamentals is ex- 
plained below. We note that the a and (3 frequencies in 
Bao.o7Ko.93Fe2As2 are 94% of those in KFe2As2, which 
is consistent with the carrier number expected from the 
composition x = 0.93. Figures 0] and \5\ show details of an- 
gular variation of dHvA oscillations in KFe2 As2 for fields 
in the (010) and the (110) plane, respectively. The an- 
gle dependences of the observed frequencies are shown in 
Fig. Efb). Figure [7] shows angular variation of dHvA os- 
cillations in Bao.o7Ko.93Fe2As2 for fields in the (010). In 
the case of Bao.o7Ko.93Fe2As2, three frequencies arising 
from materials of the pick-up coil assembly have compa- 
rable amplitudes to the frequencies from the sample: i.e., 
copper neck frequency, its second harmonic, and an un- 
known frequency, which might be due to tin from solder 
used for wiring. However, they have detectable ampli- 
tudes only for \6\ > 20° 4^ We show Fourier transforms 
of pick-up coil signals without a sample for some field di- 
rections for comparison (dotted spectra), which clearly 
indicate that the a and /3 frequencies arise from the sam- 
ple, not from the coil assembly. 



Figures [5] and [5] exemplify determination of effective 
masses for KFe2As2 and Bao.o7Ko.93Fe2As2, respectively. 
The determined masses are tabulated in Tables I and II. 
Heavy masses up to 19 m e for B \\ c are observed for 
KFe2As2. The masses in Bao.o7Ko.93Fe2As2 are nearly 
the same as the corresponding masses in KFe2As2- 



C. Fundamental frequencies 

A pair of frequencies q and e^ near FcosO = 0.3 kT 
in KFe2As2 [Figs. EJ [5j and[6^b)] can be assigned to the 
minimum and maximum orbits on a corrugated FS cylin- 
der. Their angular variation is consistent with the Ya- 
maji model as indicated in the figures. 

A complex region between F = 2 and 3 kT for B \\ c in 
KFe2As2 is shown in Fig. fTOT a) in an expanded scale. We 
identify three frequencies a/, an, and Q as fundamentals. 
All the other frequencies can be indexed with a formula 
F = F ai +nAF (n: integer), where AF = (F Ci -F Q! )/12. 
Their amplitudes are quickly suppressed as the field is 
decreased [compare the three spectra for three different 
field windows in Fig. UOf a)]. To be more quantitative, 
Fig. [11] shows the field dependence of amplitudes of var- 
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FIG. 7. (color online) (a) Angular variation of Fourier transforms (solid lines) of dHvA oscillations in Bao.o7Ko.93Fe2 As2 for the 
(010) plane. The field window is between 12 and 17.75 T for |0| < 15° and between 14 and 17.75 T for \0\ > 20°. The horizontal 
axis is FcosO. The spectra are shifted vertically so that the baseline of a spectrum measured at 9 degrees is set at 6. Three 
frequencies due to the pick-up coil (not the sample) are also observed as indicated: copper neck frequency (Cu-N), its second 
harmonic (2 x Cu-N), and an unknown frequency, which might be due to tin from solder used for wiring. These frequencies 
have detectable amplitudes only for \6\ > 20°. The spectra shown by dotted lines for some angles are Fourier transforms of 
pick-up coil signals without a sample at the same angles. The comparison indicates that the a and j3 frequencies are not from 
the pick-up coil but from the sample, (b) Fourier transforms for the field window between 12 and 14 T. Although the frequency 
resolution deteriorates because of the narrower field window, the a frequency shows up slightly clearer than in (a). 



ious frequencies for B || c in KFe2As2 in the form of 
Dingle plot. The Dingle temperatures for fundamental 
frequencies estimated from the plots are shown in Table 
I. (We can also estimate mean free paths of about 100 
and 300 nm for e and a, respectively.) It is clear that 
the non-fundamental frequencies are suppressed with de- 
creasing held more quickly than the fundamental frequen- 
cies. These observations strongly suggest that the fre- 
quencies other than ai, ah, and Q are due to magnetic 
breakdown orbits between ai and Q. We observed es- 
sentially the same spectra and held dependence in other 
samples. 

Figure [T3] shows the ai and orbits deduced from re- 
cent ARPES measurements 4£r— There are eight possible 
magnetic breakdown junctions as indicated. If area Al is 
approximately twice area A2 and hence Fax ~ 2AF and 
Fa2 ~ AF 1 , the observed magnetic breakdown frequen- 
cies can be explained: frequency n — 1 corresponds to 
F ai +F A2 ,n = 2to F at + F A1 and F ai + 2F A2 , n = 3 to 
F a , +Fai+Fa2 and F ai + 3Faz, and so on. It is interest- 
ing that the positions of the inferred magnetic breakdown 



junctions remind us of those of the octet-line nodes in the 
superconducting gap observed in a laser- ARPES study£L 

Figure UOf b) shows a frequency region near the second 
harmonic of a^h in an expanded scale. Many frequencies 
are observed above 2ai , and they (except 2ctfc) can be 
indexed with F = 2F ai + nAF (n: integer). However, a 
small peak at F = 4.40 kT can not be indexed and we 
identify it as the counterpart of Q , namely Ch ■ 

Detailed angular dependence of ai : h in KFe2As2 is 
shown in the right panels of Figs. [4] and [5j As the field is 
tilted from the c axis, the two peaks gradually approach 
and roughly merge when |0(oio)| is ~ 50° or when |#(iio)| 
is slightly larger than 40°. The amplitudes at these magic 
angles are enhanced when compared to those at neighbor- 
ing angles. These observations are basically consistent 
with the Yamaji model (dashed lines). Slight deviations 
of the peak positions from the model and extra split- 
ting of the frequencies, which is especially noticeable at 
high angles \9\ > 50°, can be explained by small error in 
and small distribution of the c-axis orientation in the 
sample. The shading in the figures illustrates how the 
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TABLE I. Experimental dHvA frequencies, effective masses, 
and Dingle temperatures in KFe2As2. 



Field direction 


Branch 


F (kT) 


m* /m e 


x' D (K) 


B || c 


Q 


0.24 


6.0(4) 


0.14(2) 




<h 


0.36 


7.2(2) 


0.18(1) 




Oil 


2.30 


6.0(2) 


0.19(2) 




ai, 


2.39 


6.5(2) 


0.18(2) 







2.89 


8.5(2) 


~0.1 




0. 


4.40 


18(2) 






P 


7.16 


19(2) 




0(oio) = 36.7° 


<-! 


0.30 


7.4(7) 






£h 


0.45 


8.4(2) 






Oil 


2.91 


7.4(2) 






o-C 
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7.8(2) 
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7.4(2) 









3.72 
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4.70 
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FIG. 8. (color online) Temperature variation of (a) ai, ah, 
and (b) /? Fourier peaks for B || c in KFe2As2- The insets 
show their amplitudes as functions of temperature. The solid 
curves are fits to the Lifshitz-Kosevich formula, from which 
effective masses are estimated. 



TABLE II. Experimental dHvA frequencies and effective 
masses in Bao.o7Ko.93Fe2As2. 



Field direction 



Branch 



F (kT) 



B 



Pi 
Ph 



2.17 
6.69 
6.77 



m* /m e 



4(1) 
18(2) 
16(3) 



m*/m. 




8.0 



FIG. 9. (color online) Same as Fig. [8] but for /3 peaks in 
Ba .o7Ko.93Fe 2 As2. 



frequencies spread if there is a ±1° error or distribution 
of the c-axis orientation. The magic angles also deviate 
from the Yamaji model and exhibit in-plane anisotropy; 
i.e, the magic angle in the (010) plane is slightly larger 
than that in the (110) plane. This indicates in-plane 
anisotropy of kp and/or more complex c-axis dispersion 
than the Yamaji model. The observed magic angles and 
their anisotropy are consistent with those observed in 
AMRO measurements^ 

The angular dependence of Cz,h [Figs. 01 [5l and[6jb)] 
indicates that the £ sheet is the most three-dimsnsional 
sheet of the Fermi surface. The maximum frequency (,h 
exhibits much more variation in F cos 9 than the mini- 
mum frequency Q. This indicates that, while the £ cylin- 
der is close to a straight one around the minimum cross 
section, it swells out locally around the maximum cross 
section. For B \\ c, the effective mass of C,h is considerably 
larger than that of Q (Table I). As the field is tilted, the 
masses of the two frequencies become nearly the same 
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FIG. 10. (color online) (a) Expanded view of Fourier trans- 
forms for B || c in KFe2As2 in a frequency range between 2 
and 3 kT. All the frequencies other than ai, olk, and C,i can be 
indexed with a formula F = F a[ + nAF (n: integer), where 
AF = (F(, — F ai )/12, as indicated by short vertical bars. A 
small peak marked by an asterisk dose not accord with the 
formula. It is however an artifact due to spectral interpola- 
tion by zero padding (the spectra shown by the solid curves 
are those interpolated by zero padding in the 1/B domain). 
The non-padded spectrum shown by dots can be consistent 
with a small n — 1 peak buried under the low-frequency tail 
of the qj, peak. Compare the three spectra for three different 
field windows: as the field is decreased, the peaks other than 
ai, ah, and Q are quickly suppressed, (b) Expanded view 
for a frequency range between 4 and 6 kT. The frequencies 
higher than the second harmonic of Q; , except for the second 
harmonic of a^, can be explained by F = 2F a[ + nAF (n: 
integer) as indicated by small vertical bars. A small peak (h 
at F = 4.40 kT can not be explained and is a fundamental 
frequency. 
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FIG. 11. (color online) Dingle plots for various frequencies 
in KFe2As2 for fl | t The frequencies indexed by n are ex- 
plained in Fig. 1101 According to the Lifshitz-Kosevich formula 
[see Eqs. (2)-(4)], this type of plot gives a straight line, the 
slope of which gives the Dingle temperature (Table I). A very 
wide field range between 3.09 and 17.65 T was used to make 
these plots. However, in order to resolve finely spaced fre- 
quencies, wide field windows are necessary for Fourier trans- 
formation, and hence a very limited number of data points is 
obtained. Note that frequencies other than a;, aj,, and Q are 
quickly suppressed as the field is decreased. 




FIG. 12. (color online) The m and Q orbits deduced from 
recent ARPES data.— ^ Eight possible magnetic breakdown 
junctions, where holes can tunnel from the m orbit to C,i (or 
vice versa) through a small energy gap, and areas AI and A2 
are indicated. 



(see the masses determined at t9(oio) = 36.7° and #(iio) 
= 40.5° in Table I). This also indicates that the defor- 
mation of the C, cylinder is local near the maximum cross 
section. 

We now turn to the /3 frequency in KFe2As2- In the 
previous paper ^ we identified this frequency with the 
third harmonic of a^ because it satisfies the relations 



Fp = 3F ah and 



"P 



3m a h 



within experimental accu- 



racy (Table I). However, we have recently observed the /3 



frequency very clearly in Bao.o7Ko.93Fe2As2 (Figs. [3] and 
[7]), where Fp is unmistakably different from 3^^. The 
effective mass of /3 in Bao.o7Ko.93Fe2As2 is close to that 
in KFe2As2 (Tables I and II). These findings have moti- 
vated us to reexamine the j3 frequency in KFe2As2, and 
we have reached the conclusion that it is a fundamental 
frequency as explained below. 

First of all, the amplitudes of the /3 frequency appear 
too big for a third harmonic (Figs. 2] and [SJ. According to 



the Lifshitz-Kosevich formula, amplitudes of harmonics 
basically decrease exponentially with the harmonic num- 
ber r. Because of the spin factor R s , r , which oscillates 
with r, the third harmonic can accidentally be compara- 
ble to or even larger than the second harmonic for some 
field directions. However, the ft frequency in KFe2As2 
has a comparable amplitude to the second harmonic of 
ah for a wide range of field directions, which seems diffi- 
cult to explain if it is the third harmonic. 

To be more quantitative, we analyze the amplitudes of 
ah, its second harmonic, and ft for B || c, using Eqs. (1)- 
(6). Since the effective mass to* and the Dingle tempera- 
ture x* D for the ah frequency have been determined from 
the temperature and field dependences of the amplitude 
(Table I), we can determine the spin-splitting parameter 
S from the amplitudes of the fundamental and second 
harmonic (Fig. HH): S = n ± 0.10294 or n ± 0.32319 
(n: integer). Using these values of 5 1 , we can estimate 
the amplitude of the third harmonic. In either case, the 
estimated amplitude is significantly smaller than the ob- 
served amplitude of the ft frequency. It thus follows that, 
although there is some contribution from the third har- 
monic of ah, the observed amplitude of ft at 9 = is 
mostly due to the fundamental frequency ft. 

There is a further evidence for the fundamental ft: the 
frequency of ft deviates from that of the third harmonic of 
ah at high angles as shown in Fig. Q3J We therefore con- 
clude that the ft frequency is a fundamental frequency. 

It is unclear at present whether the observed ft fre- 
quency is a maximum or a minimum frequency. A small 
peak at F = 7.44 kT at 8 = might perhaps be the 
counterpart (Figs. [Hand [5]). 

The experimentally determined Fermi surface in 
KFe2As2 is schematically shown in Fig. (TSl 

Lastly, we briefly discuss Bao.o7Ko.93Fe2As2. The min- 
imum and maximum frequencies of a are not resolved in 
Fig. [7J This is probably because the a oscillation is so 
weak that the quality of the signal is insufficient to resolve 
the two finely-spaced frequencies. On the other hand, the 
ft frequency splits at — (Fig. [7] or Fig|9]). If the two 
frequencies are the minimum and maximum frequencies, 
it follows that the ft cylinder is strongly two dimensional 
(AFq/Fq = 1.2%). The e and £ frequencies were not ob- 
served in Bao.o7Ko.93Fe2As2- This is probably because 
they are too weak. 

The a and ft frequencies in Bao.o7Ko.93Fe2As2 are both 
94% of those in KFe2As2 and are consistent with the 
composition x — 0.93 within experimental error. This is 
compatible with a rigid band picture. Within an effec- 
tive mass approximation, A — 2i:mi >an dEp /h . Hence, a 
dHvA frequency shift AF due to a shift in the Fermi 
energy AEp satisfies a relation AF/F — AEp/Ep. 
Therefore, if two bands have the same Ep, AF/F for 
a energy shift AEp is the same for both bands. As- 
suming the same mass renormalization for both bands, 
E F oc F/m band oc F/m*. Since F/(m*/m e ) = 0.38 kT 
for both a and ft orbits in KFe2As2, AF/F should be 
the same as observed. 
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• observed 
- fit (S = 0.10294) 
— fit (S = 0.32319) 




Harmonic number r 

FIG. 13. (color online) Amplitudes of au, its second har- 
monic, and ft for B || c in KFe2As2. The solid curves are 
fits to the first two amplitudes based on the Lifshitz-Kosevich 
formula. To illustrate the oscillatory nature of the spin fac- 
tor R a ,r, those curves are shown as if the harmonic number 
r is continuous. The observed amplitude of ft is significantly 
larger than the amplitude of the third harmonic expected from 
either fit, indicating the existence of the ft fundamental fre- 
quency. 
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FIG. 14. (color online) Fourier transform spectrum at #(oio) = 
31.7° for KFe2As2. Two additional curves are shown: they are 
the same spectrum (with the amplitude reduced by 1/10) but 
plotted against 2Fcos# and 3Fcos#, respectively, and hence 
indicate expected positions of second and third harmonics. 
The ft peak clearly deviates from the expected position of the 
third harmonic of qj. 



IV. COMPARISON TO BAND STRUCTURE 
CALCULATIONS 

The electronic band structure of KFe2As2 was cal- 
culated within the local density approximation (LDA) 
by using a full-potential linearized augmented plane- 
wave (FLAPW) method. 3 - We used the program codes 
TSPACE 50 and KANSAI-06. The experimental crystal 
structure^ including the atomic position za s of As was 
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FIG. 15. (color online) Fermi surface cross-sections observed 
via the dHvA measurements. The in-plane anisotropy is ne- 
glected. The line thickness indicates the magnitude of the 
c-axis dispersion. 



used for the calculation. Figure \W\ shows the electronic 
band structure, density of states (DOS), and FS calcu- 
lated with the spin-orbit interaction included. The DOS 
at the Fermi level {Ep) is 58.4 states/(Ry f.u.), which 
corresponds to the Sommerfeld coefficient of jband = 10.1 
mJ/K 2 mol. Note that a direct specific-heat measurement 
gives 7 = 93 mJ/K 2 mol, 52 indicating a large mass en- 
hancement of j/jband = 9 due to electronic correlations. 
As can be seen from the DOS plot, the states near Ep 
are mostly derived from the Fe 3d orbitals. Four bands 
31-34 cross Ep, and the FS consists of three concentric 
hole cylinders at the T point of the BZ (bands-32, 33, and 
34), small hole cylinders near the zone boundary (band- 
34), and a small hole pocket at Z (band-31). A previ- 
ous calculation, 5 - where Za s was relaxed by energy min- 
imization, predicted an electron cylinder at X instead of 
the hole cylinders around X. However, ARPES measure- 
ments so far have found no electron cylinder in KFe2As2 
and are qualitatively consistent with the presently calcu- 
lated FS^^£ 

For the observed dHvA frequencies, it is clear that the 
a, C, and j3 cylinders correspond to the three hole cylin- 
ders at r and the e cylinder to the small hole cylinder near 
X (Fig.[T5|). However, the quantitative agreement is poor 
as can be seen in Fig. [SJ where the calculated and exper- 
imental frequencies are shown. Calculated band masses 
ranging from 0.3 to 2.9 m e for B || c are much smaller 
than the observed effective masses (Table I), again sug- 
gesting the strong correlations. 

To see if the agreement can be improved, we have 
adjusted the band energies so that the smallest band- 
32 cylinder matches the a cylinder and that the largest 
band-34 cylinder matches the /3 cylinder. The energy of 
band-33 was adjusted to keep the total carrier number. 
The small band-31 pocket was neglected. The resultant 
frequencies are shown in Fig.^c). The large undulation 
of the £ cylinder can not be reproduced by the calcula- 
tion. 



V. DISCUSSION 

The failure of LDA calculations in KFe2As2 contrasts 
sharply with the case of BaFe2As2, where LDA calcula- 
tions give basically correct description of the Fermi sur- 
face (despite the overestimation of the antiferromagnetic 
moment)^ It seems that LDA band structure calcula- 
tions fail to predict crystal field splitting of the Fe 3d 
levels in KFe2As2- In the calculated band structure, one 
band is situated immediately above Ep at the T point 
and two nearly degenerate bands are situated above it 
[Fig. fTBT a)]. The former band is of xy character. The lat- 
ter two bands are of xz/yz character and would be truly 
degenerate if spin-orbit coupling were absent. However, 
experimental observation is that there are two similarly 
sized relatively small cylinders (a and (/) and one rela- 
tively large cylinder (/?) at the F point. This strongly 
suggests that the xz/yz bands are actually lower than 
the xy band and produce the a and £ cylinders. In fact, 
these orbital characters have been confirmed by recent 
ARPES measurements^^ 

Dynamical mean field theory (DMFT) studies have 
shown that the crystal field splitting of the Fe 3d lev- 
els may be modified if electronic correlations are treated 
beyond the LDA^i^ A recent DMFT calculation^ sug- 
gests that the order of the xy and xz/yz bands at the T 
point in KFe2As2 are inverted by electronic correlations. 
However, it should also be noted that the crystal field 
splitting is affected by hybridization between the Fe 3d 
and As 4p as well as Fe Ap states. Given the well known 
problem of band structure calculations in iron pnictides 
that structural optimization for paramagnetic states can 
not predict the correct position of As, this hybridization 
may not be calculated very accurately within the LDA, 
which would affect the crystal field splitting of the Fe 3d 
levels. This point deserves further studies. 

We can estimate the Fermi surface volume and Som- 
merfeld coefficient of the specific heat from the deter- 
mined FS within a 2D approximation (Table III). We 
have used two data sets, one at 0=0 and the other at = 
36.7°, separately for the estimation, and have obtained 
consistent results. 

The estimated FS volume corresponds to the carrier 
number of 1.01 - 1.03 holes/f.u., which is consistent with 
the stoichiometry of KFe2As2 within experimental error. 
We also compare the sizes of the observed FS cylinders 
with those obtained by ARPES^SiS and AMRO 49 mea- 
surements in Table III. The more recent ARPES data 
(Ref. |46| ) are in reasonable agreement with the present 
data. Usually AMRO measurements give quantitative 
estimates of FS sizes, but in the present case the agree- 
ment with the present data is limited, which is ascribed 
to the fact that only a small number of broad AMRO 
peaks were observed^ 

The Sommerfeld coefficient is estimated to be 82 - 94 
mJmor'K -2 , which is consistent with the direct mea- 
sure of 7 = 93 mJmoU 1 !^ 2 (Ref. |52| and also Ref. 
Wn ) and also with an estimate of 7 ~ 84 mJmol _1 K -2 
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FIG. 16. (color online) (a) Electronic band structure, (b) density of states, and (c) Fermi surface of KFe2As2 calculated with 
the spin-orbit interaction included.— Points of symmetry (solid circles) and lines of symmetry (open circles) in the Brillouin 
zone are explained in the top left figure of (c). 

TABLE III. Estimated FS volume A and Sommerfeld coefficient 7. The former is estimated from the average of a maximum 
and a minimum frequency. The latter is estimated from the average of the effective masses for the maximum and minimum 
frequencies within the 2D approximation. Two sets of data, one at B \\ c and the other at #(010) = 36.7°, are separately used 
for the estimation. FS volumes determined by ARPE S 36 ' 46 and AMRO 4 ^ measurements are also shown. 
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a There are four e cylinders in the BZ. 



in a recent ARPES study^ where effective masses have 
been determined from band dispersions near the Fermi 
level. These confirm a mass enhancement of about 9 in 
KFe2As2. Therefore, the electronic correlations are much 
stronger in KFe2As2 than in BaFe2As2, where mass en- 
hancements are 2-3. 

It is interesting to recall here that overall band-width 
renormalization factors estimated in early ARPES stud- 
ies do not vary significantly with the chemical composi- 
tion in the Bai_. r K. E Fe2As2 system: 1.5 in BaFe2As2r* 
2.7 in Ba .6Ko.4Fe 2 As2r^ and 2 in KFe 2 As 2 . 36 On the 
theoretical side, a study based on the fluctuation ex- 
change (FLEX) approximation has predicted a signifi- 
cant growth in the mass enhancement from BaFe2As2 
to KFe 2 As 2 ^ while the DMFT study of Ref. [H sug- 
gests similar mass enhancements for both BaFe2As2 and 
KFe2As 2 . These seemingly contradicting results can 
be reconciled if one notices the existence of two dis- 
tinct mechanisms contributing to the mass enhancement 
as previously pointed outi^ One is narrowing of over- 
all band widths, which can be measured as band-width 



renormalization factors in ARPES measurements and 
seems to be dealt with a DMFT fairly well. The other is 
flattening of bands near the Fermi level, which arises from 
interaction with low-energy bosonic excitations. The 
large discrepancy between the dHvA mass enhancement 
and the ARPES band-width renormalization factor in 
KFe2As2 indicates that this latter effect, which involves 
spin and orbital fluctuations, grows from BaFe2As2 to 
KF82AS2 as suggested by the FLEX study^ The exis- 
tence of strong spin fluctuations in KFe2As2 has been 
evidenced by inelastic neutron scattering^ and NMR 



measurements 
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VI. SUMMARY 

We have performed dHvA measurements on KFe2As2 
and Bao.o7Ko.93Fe2As2- For KFe2As2, we have identi- 
fied frequencies e, a, C, and j3 as fundamentals. The 
rest of observed frequencies are attributed to harmon- 
ics and magnetic breakdown orbits between ai and Q. 
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The positions of the inferred magnetic breakdown junc- 
tions between on and Q resemble those of the octet-line 
nodes of the superconducting energy gap.— With the aid 
of LDA band structure calculations, we have assigned 
the e frequency to a hole FS cylinder near the X point 
and the a, C, &ud /3 frequencies to three concentric hole 
cylinders at the Y point. On a quantitative level, how- 
ever, the agreement between the observed and calculated 
frequencies is poor. This can be attributed to the fact 
that the crystal field splitting of the Fe 3d levels are not 
correctly calculated. Effective masses are large, up to 
19 m e for B \\ c. From the dHvA data, we estimate 
the carrier concentration and Sommerfeld coefficient to 
be 1.01 - 1.03 holes/f.u. and 82 - 94 mJmol^K- 2 , re- 
spectively. They are consistent with the chemical sto- 
ichiometry and a direct measure of 93 mJmol _1 K -2 ^ 
respectively, establishing that the determined FS is com- 
plete. The large Sommerfeld coefficient, which is about 
9 times larger than a band value, indicates that KFe2As2 
is a strongly correlated metal. The discrepancy between 
the mass enhancement of 9 and the ARPES overall band- 



width renormalization of 2 (Ref. I36T ) suggests the impor- 
tance of low-energy spin and/or orbital fluctuations. For 
Bao.o7Ko.93Fe2As2, we have observed the a and /3 fre- 
quencies. The frequency change between the two com- 
pounds can be understood within a rigid-band model. 
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